Detecting photon-photon scattering in vacuum at exawatt lasers 
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In a recent paper, we have shown that the QED nonlinear corrections imply a phase correction to 
the linear evolution of crossing electromagnetic waves in vacuum. Here, we provide a more complete 
analysis, including a full numerical solution of the QED nonlinear wave equations for short-distance 
propagation in a symmetric configuration. The excellent agreement of such a solution with the result 
that we obtain using our perturbatively-motivated Variational Approach is then used to justify an 
analytical approximation that can be applied in a more general case. This allows us to find the 
most promising configuration for the search of photon-photon scattering in optics experiments. In 
particular, we show that our previous requirement of phase coherence between the two crossing 
beams can be released. We then propose a very simple experiment that can be performed at future 
exawatt laser facilities, such as ELI, by bombarding a low power laser beam with the exawatt bump. 
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I. INTRODUCTION 

Radiative corrections in Quantum Electrodynamics 
(QED) have been studied for 70 years, both theoret- 
ically and experimentally Nevertheless, in the last 
decade they have gained an increasing interest, following 
the extraordinary advancements in the fields of Quan- 
tum and Nonlinear Optics. On one hand, it has been 
noted that a fundamental uncertainty in the number 
of photons is unavoidably generated by QED radiative 
corrections [2], eventually competing with the experimen- 
tal errors. On the other hand, the interchange of virtual 
electron-positron pairs can produce classically forbidden 
processes such as photon-photon scattering in vacuum. 

Although it is a firm prediction of QED, photon- 
photon scattering in vacuum has not yet been detected, 
not even indirectly. The rush for its discovery is then 
wide open. A hammer strategy will be to build a photon- 
photon collider (3}, based on an electron laser producing 
two beams of photons in the MeV range (i.e., having 
wavelengths in the range of a few fm). This will maxi- 
mize the cross section for the process. A second approach 
will be to perform experiments using ultrahigh power op- 
tical lasers, such as those that will be available in the 
near future^], in such a way that the high density of 
photons will compensate the smallness of the cross sec- 
tion. In this case, the photon energies are well below the 
electron rest energy, and the effect of photon-photon col- 
lisions due to the interchange of virtual electron-positron 
pairs can be expressed in terms of the effective Euler- 
Heisenberg nonlinear Lagrangian^, Q. This modifies 
Maxwell's equations for the average values of the elec- 
tromagnetic quantum fields [6( and affects the properties 
of the QED vacuumQ. 

Ultra intense photon sources are available thanks to 
the discovery of chirped pulse amplification (CPA)[9( in 
the late 80's and optical parametric chirped pulse am- 
plification (OPCPA) in the 90's. These techniques 



opened the door to a field of research in the bound- 
ary between optics and experimental high-energy physics, 
where lots of novelties are expected to come in the next 
years. In fact, several recent works propose different con- 
figurations that can be used to test the nonlinear opti- 
cal response of the vacuum, e.g. using harmonic gen- 
eration in an inhomogencous magnetic field 11], QED 
four- wave mixing (T^|. resonant interactions in microwave 
cavities [l~3|, or QED vacuum birefringence [T3 which can 
be probed by x-ray pulses[lj|, among others[l6|. 

In a recent Letter we have shown that photon- 
photon scattering in vacuum implies a phase correction 
to the linear evolution of crossing electromagnetic waves, 
and we have suggested an experiment for measuring this 
effect in projected high-power laser facilities like the Eu- 
ropean Extreme Light Infrastructure (ELI) project [T^] 
for near IR radiation. Here, we provide a more complete 
analysis, including a full numerical solution of the QED 
nonlinear wave equations for short-distance propagation 
in a symmetric configuration. The excellent agreement 
of such a solution with the result that we obtain us- 
ing a perturbatively-motivated Variational Approach is 
then used to justify an approximation that can be ap- 
plied in a more general case. This allows us to find the 
most promising configuration for the search of photon- 
photon scattering in optics experiments. In particular, 
we show that our previous requirement of phase coher- 
ence between the two crossing beams can be released. 
We then propose a very simple experiment that can be 
performed at future exawatt laser facilities, such as ELI, 
by bombarding a low power laser beam with the exawatt 
bump. The effect of photon-photon scattering will be de- 
tected by measuring the phase shift of the low power laser 
beam, e.g. by comparing it with a third low power laser 
beam. This configuration is simpler, and significantly 
more sensitive, than the one proposed in our previous 
paper. Even in the first step of ELI, we find that the 
resulting phase shift will be at least A<I> « 2 x 10~ 7 rad, 
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which can be easily measured with present technology. 
Finally, we discuss how the experimental parameters can 
be adjusted to further improve the sensitivity. 

The paper is organized as follows: 

In Section II, following Ref. [13], we present the non- 
linear equations that replace the linear wave equation 
when the QED vacuum effect is taken into account. 

In Section III, we consider a linearly polarized wave, 
describing the scattering of two counterpropagating plane 
waves that have the same intensity and initial phase at 
the origin (symmetric configuration). In this case, we 
perform a numerical simulation for short-distance prop- 
agation, which provides the first known solution to the 
full QED nonlinear wave equation. 

In Section IV, we study the same configuration as in 
Sect. Ill by using a Variational Approximation. This 
allows us to obtain an analytical solution which is valid 
even for a long evolution. After showing that this varia- 
tional solution is in good agreement with the numerical 
simulation of Sect. Ill, we give a method to study differ- 
ent configurations that cannot be integrated numerically 
in a simple and direct way. 

In Section V, we apply our Variational Approximation 
to the case of the scattering of a relatively low power 
wave with a counterpropagating high power wave, both 
traveling along the z axis. We allow an arbitrary initial 
phase relation between the two waves, and obtain an an- 
alytical solution showing a phase shift of the low power 
wave. 

In Section VI, we discuss the possibility of detect- 
ing photon scattering by measuring the phase shift of 
two counterpropagating waves at ELI. We show that an 
asymmetric configuration, in which only one high power 
laser beam is used and the phase is measured on the low 
power laser beam, is not only simpler to realize, but it 
also gives a better sensitivity than the symmetric con- 
figuration. The resulting proposed experiment will then 
allow to detect photon-photon scattering as originated 
by QED theory. 

In Section VII, we resume our conclusions, and dis- 
cuss why we think that our proposed experiment will be 
the simplest and most promising way to detect photon- 
photon scattering with optical measurements in the near 
future. 



II. THE NONLINEAR EQUATION FOR 
LINEARLY POLARIZED WAVES IN VACUUM 

In this section we introduce the equations that govern 
the evolution of the electromagnetic fields E and B when 
the QED effects are taken into account. We assume that 
the photon energy is well below the threshold for the pro- 
duction of electron-positron pairs, 2m e c 2 ~ lMeV. This 
means that we will only consider radiation of wavelengths 
A > 2 x lCP 13 m, which is always the case in optical ex- 
periments. In this case, the QED effects can be described 
by the Euler-Heisemberg effective Lagrangian density 



C = C + CC Q = Co + £ 
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Co = ^ ( E 2 - c 2 B 2 
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(1) 



(2) 



the linear Lagrangian density and eo and c the dielectric 
constant and the speed of light in vacuum, respectively. 
As it can be appreciated in Eq. (Q]) , QED corrections are 
introduced by the parameter 
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6 . 7 x 10" 30 ^-. 



(3) 



This quantity has dimensions of the inverse of an en- 
ergy density. This means that significant changes with 
respect to linear propagation can be expected for val- 
ues around |££o| ~ 1, corresponding to beam fluxes with 
electromagnetic energy densities given by the time-time 
component of the energy-momentum tensor 



dC 



[ oo 



d(d t A) 



d t A - C > 2/£ ~ 3 x 10 2y J/m 
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While such intensities may have an astrophysical or cos- 
mological importance, they are not achievable in the lab- 
oratory. The best high-power lasers that are being pro- 
jected for the next few decades will be several orders of 
magnitude weaker, eventually reaching energy densities 
of the order p ~ 10 23 J/77i 3 [l]. Therefore, we will study 
here the "perturbative" regime, in which the non-linear 
correction is very small, |££q| -C 1. As we shall see, even 
in this case measurable effects can be accumulated in the 
phase of beams of wavelength A traveling over a distance 
of the order A|££o| -1 - Thus, current sensitive techniques 
could be used to detect traces of QED vacuum nonlin- 
earities. 

Once the electromagnetic fields are expressed in terms 
of the four-component gauge field = (A , A) as B = 
V A A and E = -cV^4° - the equations of motion 
are given by the Variational Principle: 



ST 



= 0, 



(•5) 



where T = J Cd 4 x is the QED effective action. Instead of 
studying the resulting equations for the fields E and B, 
that can be found in the literature [H 0], for the present 
purposes it is more convenient to consider the equations 
for the gauge field components A^ . In general, these four 
equations cannot be disentangled. However, after some 
straightforward algebra it can be seen that they admit 
solutions in the form of linearly polarized waves, e.g. in 
the x direction, with A = and A = (A, 0, 0), provided 
that: i) the field A does not depend on the variable x 
(a transversality condition) and it) A(t, y, z) satisfies the 
single equation: 
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+ £e c 2 {d^Ad v Ad v A + 2d ll Ad l tAd v A) = 0, 

where we have used the convention g^ v = 
diag(l, — 1, — 1, — 1) for the metric tensor. In non- 
relativistic notation, Eq. © can be written as 

= d 2 y A + d 2 z A - d 2 t A + £e c 2 { (7) 
[ (d t A) 2 - 3(d y A) 2 - (d z A) 2 ] d 2 A+ 
[ (d t A) 2 - {d y Af - 3(d z A) 2 ] d 2 A- 
[ 3(d t A) 2 - (d y A) 2 - (d z A) 2 ] d 2 A+ 
4 (d z Ad t Ad z d t A - d z Ad y Ad z d y A + d y Ad t Ad y d t A) } , 

where d y = ^, 3 Z = £ and d t = Jj. 

Hereafter, we will restrict our discussion to this case 
of linearly-polarized solution. The orthogonality relation 
E • B = is then automatically satisfied, and the effec- 
tive Lagrangian Eq. ([I]) reduces to C = £o(l + ££o). 
Note that the plane- wave solutions of the linear Maxwell 
equations, such as A cos (k • r — ut), where A is a con- 
stant, k = (0,k y ,k z ) and uu = c|k|, are still solutions of 
Eq. ([5]). However, we expect that the non-linear terms 
proportional to £, due to the QED correction, will spoil 
the superposition principle. 

III. A NUMERICAL SOLUTION OF THE FULL 
NONLINEAR EQUATION (SYMMETRIC 
CONFIGURATION) 

In general, the numerical solution of a nonlinear wave 
equation such as Eq. |7|) is a formidable problem. In 
fact, it is of an higher order (and is much more compli- 
cated) than the nonlinear Shrodinger equation in three 
dimensions, which is nevertheless a highly nontrivial and 
rich system[20] . There are two main difficulties in dealing 
with Eq. ([7]). First of all, the direct numerical integration 
is practically impossible when the integrating interval is 
much larger than the wavelength. Second, one has to find 
convenient boundary conditions. 

In this section, we will study a particular configura- 
tion for which a numerical solution of the full nonlinear 
equation Eq. ((7J) can be found. Let us first consider 
two counter-propagating plane waves that travel along 
the z-axis, for simplicity having the same phase at the 
space-time origin. The corresponding analytical solution 
of the linear wave equation (that can be obtained by set- 
ting £ = in Eq. (JTJ) ) would be 

A 

j4ij n (t, z) = — [cos(kz — uut) + cos(fcz + ut)] (8) 
~ Acos(cut) cos(fcz), 

where A is a constant amplitude, k = (0, 0, fc) is the 
wave vector, and u> = ck is the angular frequency. It is 
easy to see that Eq. ([5]) can also be considered as the 



analytical solution of the linear wave equation satisfying 
the boundary conditions 

A(t,0) = Acoa(ujt), (9) 

O z A(t,0) = 0, 

where for convenience we chose as the integration inter- 
val a 'small' cuboid of time-dimension 2h/oj and space- 
dimension 2-7r/fc. With such a choice, the linear wave 
equation can also be integrated numerically, obtaining a 
result that coincides (within the numerical error) with 
the analytical one. 

Of course, in the case of the linear equation this result 
is trivial. However, we will use it here as a guide in 
order to find a solution of the nonlinear Eq. ([7]). In fact, 
although we do not know any analytical solution of the 
latter equation, we can find its numerical solution that 
satisfies the same boundary conditions of Eq. © using 
the same small cuboid as the integration interval. 




FIG. 1: Numerical solution A nnm /A of Eqs. (O and (J9j for 
S;£oA 2 u> 2 = 0.01, as a function of the adimensional time and 
space coordinates r = cut and f = kz. 

With these conditions, the result of the numerical in- 
tegration of Eq. ([7]) is shown in Figs. [1] and [3 for a 
choice of parameters such that £,eoA 2 oj 2 — 0.01. The 
corresponding energy density, obtained as 

P = Too (10) 
= ^(E 2 + c 2 B 2 ) + |e 2 (E 2 - c 2 B 2 ) (3E 2 + c 2 B 2 ) , 

is plotted in Fig. [3l 

As far as we know, this is the first time that a numer- 
ical solution of the full nonlinear wave equation Eq. (J7J) 
is obtained. Of course, it corresponds to a very simple 
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FIG. 2: Zero-time plot of A nunl /A, as a function of the adi- 
mensional space coordinate £ = kz. 




FIG. 3: Energy density corresponding to the solution A nuln , 
as a function of the adimensional time and space coordinates 
r = ujt and C, = kz. 



particular case and a short space-time evolution. How- 
ever, we will see later that this solution will allow us to 
reach very important conclusions and prove the viability 
of useful approximation methods. 

For the moment, we note that our numerical solution, 
as given in Figs. [2] and O shows an oscillatory be- 
haviour. However, the time-average of the energy den- 
sity turns out to be approximately constant along the 
^-evolution, giving £p = ^£ f^^f u pdt = 0.00252. As we 
have discusses above, this value of the product £/? is sev- 
eral order of magnitude larger than what can be achieved 
in the laboratory in the next decades. Of course, we will 
use realistic values of p when we will present our propos- 
als of experiments in the last sections. For the moment, 
it is interesting to note that even such an enormous en- 
ergy density is still small enough so that the effect of 
the nonlinear terms gives a small correction to the lin- 



ear evolution in the short distance. In fact, we see from 
Fig. 2] that the relative difference between the linear and 
the nonlinear evolution, in our integration interval which 
is of the order of the wavelength, is of the order of few 
percent, i.e. of the same order than the adimensional pa- 
rameter £eo-4 2 ^ 2 = 0.01. This result is not surprising, 
and will provide a justification for the perturbatively- 
motivated variational approach that we will use in the 
next sections. 




FIG. 4: Relative error (A num — An n ) /A of the linear approx- 
imation, as a function of the adimensional time and space 
coordinates r = uit and f = kz 



However, from Fig. Hlitself, we can also appreciate that 
the difference between the linear and nonlinear behaviour 
tends to increase along the z evolution, so that it can 
be expected that it will eventually become large after a 
distance much larger than the wavelength. In the next 
section, we will see an analytical argument that confirms 
this expectation. 
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FIG. 5: Detail of the zero-time functions A nuln / A (upper 
curve) and An n /A (lower curve), for values of the adimen- 
sional space coordinate £ = kz close to the second zero. 
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Finally, in Fig. [5] we compare the z-evolution of 
the solutions ^4 num and An n in a greater detail for val- 
ues of z around the second zero of the solutions. We 
see that A num anticipates A\ ln , and the corresponding 
phase shift can be evaluated numerically if we define an 
effective wave vector component k z by computing the 
value zq corresponding to ^4ii n (0,zo) = 0, and setting 
k z zo — 3tt/2. The numerical determination of the zero 
gives Co = kz = 4.68885, so that k z = 1.0050fc. We will 
provide a full explanation for this result in the following 
section. 



IV. VARIATIONAL APPROXIMATION 
(SYMMETRIC CONFIGURATION) 

Although it can be considered as an interesting 
achievement due to its simplicity and lack of previous ap- 
proximations, in practice the numerical solution that has 
been discussed above can only be obtained in the spe- 
cial configuration of two counter-propagating, in-phase 
waves, and for short propagation (of the order of the 
wavelength). In this section, we will consider the same 
configuration, and we will look for a variational approx- 
imation that will provide an analytical solution which is 
valid even for a long evolution. After showing that this 
variational solution is consistent with the numerical one 
of the previous section, we will obtain a method to study 
different configurations that cannot be integrated numer- 
ically in a simple and direct way. 

We will thus consider again the two counterpropagat- 
ing waves, that would be described by Eq. {SJ in the 
linear case. Note that any of the two crossing waves, 
Y cos(fcz — tot) and =j cos(fcz + wt), when taken alone, 
would be a solution of both the linear and non-linear 
equations, provided that u> = ck. However, their super- 
position would only solve the linear equations of motion. 

As we discussed in Ref. [l7|, in all the experimental 
configurations that can be studied in the present and 
the near future, the product £/? will be so small, that 
the non-linear correction will act in a perturbative way, 
progressively modifying the form of A(t, y, z) as the wave 
proceeds along the z direction. We will therefore make 
the ansatz 

A = A [a{z) cos(fcz) + f3(z) sin(A;z)] cos(wi), (11) 




The expression that is obtained by this procedure is 
still quite complicated, due to the presence of the trigono- 
metric functions of multiples of kz. However, as we shall 
see below, when £p <C 1 this z-dependence is much faster 
than that of the envelop functions a and (3. Therefore, it 
is a very good approximation to perform the z integral in 
two steps: first, we average over a period, treating a and 
(3 as if they were constant. In this way, we get rid of the 
trigonometric functions. At this point, we allow again 
the z dependence of the envelop functions. If we call T 
the average action that is obtained by this procedure, we 
will then minimize it with respect to the functions a and 
/3 by solving the equations <5r /5a = y ST /6[3 = 0. 

After a long but straightforward computation, and ne- 
glecting the non-linear terms involving derivatives of the 
functions a and /3, since they give a smaller contribution 
as we discussed in Ref. [l7| . these two equations can be 
written as 



— +/3' + x(a 2 + /? 2 )a = 0, (13) 

and 



^-«' + X (« 2 +/5 2 )/?-0, (14) 

where x = £,eoA 2 c 2 k 3 /2 is a parameter that describes the 
leading nonlinear effects. 

Now, in the perturbative regime in which \ is small, 
we can also neglect the second derivatives in the previous 
equations so that an (approximate) analytical solu- 
tion for Eqs. (fT3|) and (fT4|) is given by a(z) = cos(xz) 
and (3(z) — — sin(xz). Substituting this result into the 
variational ansatz (|11[) . and using elementary trigonom- 
etry, we obtain 

A = Aco$(wt) cos[(fc + x)z], (15) 

As a result, we find that the phase of the wave is shifted 
by a term 



allowing for the generation of the other linearly- 
independent function sin(fcz) (we will take a(0) = 1 and 
(3(0) = 0). Note that the invariance of Eq. ([7]) under time 
reversal guarantees that, if the initial behavior is propor- 
tional to cos(a;i), which is even under time inversion, no 
uneven term proportional to sin(wi) will be generated. 

According to the Variational Method, the best choice 
for the functions a(z) and (3(z) corresponds to a local 
minimum of the effective action T, after averaging out 
the time-dependence as follows: 



Acj) = xz. (16) 
Note that the solution of Eq. (| 1 5[) can also be written 

as 



A = - cos[(fe + x)z - wt] + - cos[(fc + x )z + tot], (17) 

therefore we see that each of two scattering waves is phase 
shifted according to Eq. (HHJ), due to the crossing with 
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the other wave. A similar behavior can also be found in 
a nonlinear medium [2l|. although the analogy cannot be 
pushed too far, as we have discussed in Ref. [l7(, where 
we have shown that the vacuum does not present the 
usual AC Kerr effect. 

Note also that the result of Eq. (fT"5|) can be stated 
equivalently by defining a wave vector as k z — k + x, that 
satisfies a modified dispersion relation, uj — c{k z — x) ■ 

In order to make definite numerical predictions, it is 
convenient to express the amplitude A in terms of the 
energy density, as given by Eq. (flU)) . In our pertur- 
bative regime, after time and space average, this gives 
p ~ eoA 2 uj 2 /4: with a very good approximation. Thus we 
get x — 2£p&, so that the phase shift accumulated after 
a distance Az is 

Acj) ~ 2£pkAz (18) 

We see now that the slow varying envelop approxima- 
tion was justified as far as x ^ k. 

For the same choice of parameters that was considered 
in the previous section, £p ~ £eoA 2 uj 2 /4: = 0.0025, in 
agreement with the numerical simulation. Thus we get 
X = 0.005fc, which is two order of magnitude smaller 
than k. Our approximations can then be expected to be 
reasonably good even for the extremely large of p that we 
have chosen here. Note also that this corresponds to a 
value k z = k + x — l-005fc, in agreement with the result 
that we obtained from the numerical simulation in the 
previous section. 




FIG. 6: Relative error (A nU m — Avar) I A of the variational 
approximation, as a function of the adimensional time and 
space coordinates t = cut and ( = kz. 

In Fig. [5] we compare the corresponding analytical 
solution A var , given in Eq. (|15p , with the numerical solu- 
tion A num of the QED wave equations that we have found 
in the previous section. Comparing with Fig. 2J we see 
that the variational solution is an order of magnitude 



closer to the numerical simulation than the linear evolu- 
tion, Eq. ([5]) , that was obtained by completely neglecting 
the nonlinear terms. This is a significant improvement in 
such a short distance. However, according to the previ- 
ous discussion, in the perturbative regime corresponding 
to small values of the product £p, the variational solu- 
tion is expected to be a good approximation even when 
a longer propagation distance is considered along the z- 
axis. In fact, this expectation is confirmed by Fig. [6] 
itself, that shows that the error of the variational solu- 
tion oscillates without any substantial increment in the 
integration interval. As we have observed in the previ- 
ous section, this was not the case for the linear solution, 
which was incresing its error even in the short distances. 
Of course, this is due to the fact that it does take into 
account the phase shift, which is the leading effect due to 
the nonlinear terms according to our Variational Method. 
On the other hand, the agreement of the variational so- 
lution with the numerical simulation can be used as an 
additional, a posteriori justification for our analytical ap- 
proach. 

V. VARIATIONAL APPROXIMATION 
(ASYMMETRIC CONFIGURATION) 

In the previous section, we have proved the reliabil- 
ity of a perturbatively-motivated Variational approach 
in the search for approximate solutions of the QED non- 
linear wave equation, Eq. ([7]) In that case, we studied a 
symmetric configuration in order to compare the result 
with the numerical integration of the full equation. With 
this strong justification in mind, we can now apply the 
Variational Method to different configurations, that do 
not allow for a direct integration of the full Eq. ([7]). 

In this section, we will study the case of the scattering 
of a relatively low power wave with a counterpropagating 
high power wave, both traveling along the z axis. In this 
case, the main effect will be the modification of the low 
power beam due to the crossing with the high power one, 
which would be unaffected in a first approximation. 

Using our variational approach, we will describe the 
high power field by a wave polarized in the x direction, 
having the x-component of the vector potential given by 

A h = Acos(kz + Lut + ip), (19) 

where (p is an arbitrary phase, describing the unknown 
phase difference between the two counterpropagating 
beams. Such a beam is shot against a low power wave 
that, in the absence of the nonlinear QED terms, would 
be described by a ^-component of the vector potential 
given by Ai(t, z) — acos(kz — ujt), where the constant a 
is related to the intensity of the beam. In other words, 
the only condition is that the two waves have the same 
frequency. Photon-photon scattering will then generate 
a z-dependence of a and an additional term proportional 
to sin(fcz — ujt) in the low power wave, so that 
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Ai(t, z) = a(z) cos(kz - tut) + /3(z) sin(kz - cut), (20) 

with the initial condition a(0) = ao and (3(0) = 0. Ne- 
glecting the effect of the low power beam on the high 
power wave, we are then lead to the ansatz 



A{t,z) = A h (t,z)+A l {t,z). (21) 

According to the Variational Method, we require that 
the functions a(z) and (3(z) correspond to a local min- 
imum of the effective action T, after averaging out the 
time-dependence as described in Eq. (fT2|) . After av- 
eraging over the fast dependence on z, as discussed in 
the previous section, and neglecting again the non-linear 
terms involving derivatives of the functions a and (3, we 
get the following equations 



— +P' + m = 0, (22) 

and 



^ - d + 77/3 = 0, (23) 

where ry = 2£eoA 2 c 2 k 3 . Note that Eqs. ((22j) and (23|) 
are linear, and they do not show any dependence of the 
initial phase difference ip between the two beams. Ne- 
glecting the second derivatives, for the same reason that 
were discussed in the previous section, we get the solu- 
tion a(z) = ag cos(r]z) and /3(z) = — aosin(^z). Note 
also that this solution is valid for any value of the am- 
plitude ao of the low power wave, provided that it is 
much smaller than the amplitude of the high power beam, 

Kl < 1-41- 

As a result, we find that the low power beam becomes 



Ai(t,z) = aocos[(k + r))z- ut], (24) 

Taking into account that the average energy density of 
the high power wave is p = eo»4 2 <x> 2 /2, we can compute 
the phase shift accumulated by the low power wave after 
a distance Az as 



A$ = r]Az ~ 4£pfcAz. (25) 

We stress that this result only depends on the energy 
density of the high power wave, as far as it much larger 
than that of the low power wave. The initial phase dif- 
ference is found to be irrelevant. 



VI. PROPOSAL OF AN EXPERIMENT 

We can now discuss the possibility to test the non- 
linear properties of the QED vacuum by measuring small 
phase changes in one of two crossing laser beams at a very 
high power laser facility. In the previous sections, we have 
studied the leading QED nonlinear effect on two counter- 
propagating waves traveling in the z-direction. We have 
considered two different configurations: a symmetric one, 
in which both waves are high power waves with the same 
initial phase and amplitude; and an asymmetric config- 
uration, in which only one the two waves is high power, 
and they do not need to be in phase. Of course, the 
second configuration is easier to be produced experimen- 
tally. Moreover, as we shall discuss below, it leads to 
a phase shift which is (at least) twice higher than that 
which could be obtained from the symmetric configura- 
tion, for the same experimental facility. 

In fact, the total energy density achievable in the sym- 
metric configuration is the space-time average due to 
both waves, that have to be obtained from the same origi- 
nal high power pulse through a beam splitter. Therefore, 
even if we neglect the energy loss when the beams split, 
the value of p in Eq. ([23)) is roughly the same than that 
of the single high power beam of Eq. (fT8"j). Therefore, the 
asymmetric configuration of section V produces a phase 
shift which is roughly twice than that of the symmetric 
configuration that we have discussed. 

Therefore, we will propose the following experimental 
setup, which in our opinion will provide the simplest and 
most effective way to look for optical effects of photon- 
photon scattering in future exawatt laser facilities. A 
common laser pulse is divided in two beams, A and B, 
one of which (say A) crosses at a 180° a very high power 
beam. As a result, the central part of the distribution of 
the beam A has acquired a phase shift A<& with respect 
to beam B, that has propagated freely. In a experiment 
corresponding to the parameters of the ELI project in its 
first step we have pulses of wavelength A = 800nm, inten- 
sity / = 10 29 Wm~ 2 and duration At = 10/s which are 
focused in a spot of diameter d w 10p.m. From Eq. (|25[) . 
this results in a phase shift A<E> « 2 x 10~ 7 rad for beam A, 
which can be resolved comparing with the beam B which 
was not exposed to the effects of QED vacuum. Cur- 
rent techniques like spectrally resolved two-beam cou- 
pling, which can be applied for ultrashort pulses [13], can 
be used to this purpose. 

It is interesting to note that the sensitivity of this 
method for the detection of photon-photon scattering 
may be enhanced by a suitable choice of the combina- 
tion of the intensity /, the wavelength A and the time 
duration At that enter in Eq. (|25[) . In fact, taking into 
account that Az ~ cAt, Eq. (|23|) implies that the most 
sensitive experimental configuration will be that having 
the maximum value of the combination I At/ A. 

Comparing to other alternatives like x-ray probing of 
QED birefringence, our system does not need an extra 
free electron laser and the power requirements of the sys- 
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tern are only one order of magnitude higher. Moreover, 
the measurement of the ellipticity and the polarization 
rotation angle in birefringence experiments is not yet pos- 
sible with current technology. Other techniques like four- 
wave mixing processes (12j require the crossing of at least 
three beams, with the corresponding alignment problems 
and the rest of the requirements are similar to our pro- 
posal. Moreover, the present result also improves signif- 
icantly the one arising for two high-power waves that we 
discussed in Ref. fi"7l ]. 

VII. CONCLUSIONS 

In this paper, we have studied the nonlinear wave equa- 
tions that describe the electromagnetic field in the vac- 
uum taking into account the QED corrections. In par- 
ticular, we have studied the scattering of two waves in a 
symmetric and an asymmetric configuration. In the first 
case, we have found a full numerical solution of the QED 
nonlinear wave equations for short-distance propagation. 
We have then perfomed a variational analysis and found 
an analytical approximation which is in very good agree- 
ment with our numerical solution, but can also be used 
for long distance propagation. We have then studied an 
aymmmetric configuration corresponding to the head-on 
scattering of a ultrahigh power with a low power laser 
beam and argued that it is the most promising configu- 



ration for the search of photon-photon scattering in op- 
tics experiments. In particular, we have shown that our 
previous requirement of phase coherence between the two 
crossing beams can be released. We have then proposed 
a very simple experiment that can be performed at fu- 
ture exawatt laser facilities, such as ELI, by bombarding 
a low power laser beam with the exawatt bump. Photon- 
photon scattering will then be observed by measuring the 
phase shift of the low power laser beam, e.g. by compar- 
ing with a third low power laser beam. This configura- 
tion is simpler, and significantly more sensitive, than that 
proposed in our previous paper. Even in the first step of 
ELI, we have found that the resulting phase shift will be 
at least A$ w 2 x 10~ 7 rad, which can be easily measured 
with present technology. Finally, we have discussed how 
the experimental parameters can be adjusted to further 
improve the sensitivity. 
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